Noise suppression due to long-range Coulomb interaction: 
Crossover between diffusive and ballistic transport regimes 
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( |cond-mat/9703123| , Semicond. Science and Techno!., v. 12, 1053 (1997)) 

We present a Monte Carlo analysis of shot-noise suppression due to long-range Coulomb interaction 
in semiconductor samples under a crossover between diffusive and ballistic transport regimes. By 
varying the mean time between collisions we find that the strong suppression observed under the 
ballistic regime persists under quasi-ballistic conditions, before being washed out when a complete 
diffusive regime is reached. 



It has been recognized that shot noise and thermal 
noise are not different sources of noise, but just special 
limits of a more general concept (see discussion by Lan- 
dauer Qj), which can be treated in a unified way [jl],||. 
This point has been clarified by studying noise properties 
of mesoscopic systems such as quantum point contacts, 
quantum wires, quantum dots, etc || in which both bal- 
listic (X p 3> L) and diffusive (X p <C L) regimes of carrier 
transport are accessible (here X p is the mean free path 
and L is the device length). 

A matter of particular interest is the suppression of 
shot noise resulting from the correlation among carriers 
in their transmission through a device (see some recent 
experiments @,||,^,@,D and the status of theoretical devel- 
opments ||). Among the different mechanisms of noise 
suppression most attention has been paid to the corre- 
lations imposed by the Pauli principle (Fermi statistics) 
under degenerate conditions. In this case, the reduction 
of noise has been found for both ballistic || and diffusive 
regimes [fl0|], and the crossover between them has been 
studied theoretically by de Jong and Beenakker [jl]| . 

However, as emphasized by Landauer [JjJ, apart from 
Pauli correlations, Coulomb interaction between carri- 
ers may also be a source of noise suppression. By us- 
ing an ensemble Monte Carlo simulator self-consistently 
coupled with a Poisson solver (PS) we have calculated 
the shot-noise suppression associated with long-range 
Coulomb interaction in the case of perfect ballistic trans- 
port . Here we address the case of low-density carrier 
concentration typical of semiconductors, thus avoiding 
any effect due to statistical degeneracy. Under this con- 
dition we have found that this suppression is stronger 
as space-charge effects become more important, and can 
achieve higher levels for increasing values of the param- 
eter A = L/Ld c , where Ljjc = \/ £ksT / 'q 2 n c is the De- 
bye length corresponding to the carrier concentration at 
the injecting contact (T is the lattice temperature, ks 
is Boltzmann constant, e is the dielectric constant, q is 
the electron charge, and n c is the carrier concentration 



at the contact). Therefore, the longer the device length, 
the larger the shot-noise suppression is expected to be, 
provided that transport remains ballistic. However, with 
increasing the device length (or lattice temperature) car- 
rier transport actually becomes diffusive and the disap- 
pearance of noise suppression associated with Coulomb 
correlation is expected. Note that in the case of met- 
als all charge fluctuations are screened out and Coulomb 
correlations play no role in the dc transport regime. It is 
the aim of the present letter to investigate the behavior of 
noise suppression due to the long-range Coulomb interac- 
tion among carriers in the crossover region from ballistic 
to diffusive transport regimes in nondegenerate semicon- 
ductors. We will show that noise suppression still remains 
important under quasi-ballistic conditions, before being 
washed out by the action of scattering mechanisms. 

For the calculations we consider the following simple 
model: a lightly doped active region of a semiconductor 
sample sandwiched between two heavily doped contacts 
(of the same semiconductor) injecting carriers into the 
active region. The contacts are considered to be ohmic 
(the voltage drop inside them is negligible) and they 
remain always in thermal equilibrium. Thus electrons 
are emitted from the contacts according to a thermal- 
equilibrium Maxwell-Boltzmann distribution and they 
move inside the active region, following the scmiclassical 
equations of motion, by undergoing elastic and isotropic 
scattering. The fluctuating injection rate at the con- 
tacts, which is associated with the diffusion current of a 
homojunction, is taken to follow a Poissonian statistics. 
Accordingly, the time between two consecutive electron 
injections is generated with a probability per unit time 
P(t) = re _rt , where T = ^n c vthS is the injection rate, 
Uth = ^/SksT/nra is the thermal velocity, S is the cross- 
sectional area of the device and m is the electron effective 
mass. The simulation is one-dimensional in real space 
(Poisson's equation is solved only in one dimension) and 
three-dimensional in momentum space. 

For the simulations we have used the following set of 



parameters: T — 300 K, m — 0.25 too (mo being the 
free electron mass), relative dielectric constant £=11.7, 
L = 2000 A and n c — 4 x 10 17 cm -3 , much higher than 
the sample doping (here taken of 10 11 cm -3 , but the same 
results are obtained up to 10 15 cm~ 3 ). The above set of 
values yields for the dimensionless parameter A = L/Ld c , 
which characterizes the importance of the electrostatic 
screening, the value A = 30.9. Along with A, we intro- 
duce the dimensionless ballistic parameter £ — vtwr/L 
which characterizes the crossover between diffusive and 
ballistic transport regimes. The average time between 
collisions in the bulk r is assumed to be independent of 
energy [Q, and it is varied from 10~ 15 s to 10" 11 s, so 
that both regimes of carrier transport, ballistic (f > 1) 
and diffusive (I « 1), are covered. For the PS, a time 
step of 2 fs and 100 meshes in real space are taken. 

We apply a dc voltage and calculate the time-averaged 
current / and the current autocorrelation function Ci(t) 
by means of the ensemble Monte Carlo simulator self- 
consistent ly coupled with the PS 0. We stress that in 
our approach the number of electrons N inside the sam- 
ple is not required to be fixed. Carriers are injected at 
x = and x = L inside the active region of the device 
according to P{t). When a carrier exits through any of 
the contacts it is cancelled from the simulation statistics, 
which accounts only for the carriers which are inside the 
device at the given time t. Thus, N(t) fluctuates in time 
due to the random injection from the contacts and we can 
evaluate both the time-averaged value (N) and its fluctu- 
ations 8N(t) — N(t) — (N). To illustrate the importance 
of the effects associated with the long-range Coulomb 
interaction, we provide the results for two different sim- 
ulation schemes. The first one involves a dynamic PS, so 
that any fluctuation of space charge arising because of the 
random injection from the contacts causes a redistribu- 
tion of the potential which is self-consistently updated by 
solving the Poisson equation at each time step during the 
simulation. The second scheme makes use of a static PS, 
so that only the stationary potential profile is calculated; 
once the steady state is reached, the PS is switched off 
and carriers move in the frozen non-fluctuating electric 
field profile. Both schemes give exactly the same average 
current and steady-state spatial distributions of all the 
quantities, but the noise characteristics are quite differ- 
ent. 

By varying the ballistic parameter £ through r we find 
that the steady-state spatial profile of the potential does 
not change significantly. It always exhibits a minimum 
near the cathode due to the presence of space charge 
whenever the current does not saturate. This minimum 
acts as a potential barrier for the electrons moving be- 
tween the contacts, so that those electrons which do not 
have enough energy to pass over the barrier are reflected 
back to the contacts. The most important fact is that 
the barrier height and, as a consequence, the transmis- 
sion through it depends on the current, which is crucial 
in calculating the noise characteristics. 

Fig. [l] presents the results for the low-frequency spec- 
tral density of current fluctuations Si = 2 Ci(t)dt as 
a function of I for an applied voltage of U — 40ItBT/q. 



For the present value of A this voltage corresponds to 
the maximum suppression of noise, reaching a value 
of 0.045. The results are normalized to 2ql s , where 
I s = qT = hqn c VthS is the saturation current (maximum 
current that our contacts in thermal equilibrium can pro- 
vide). Both static and dynamic PS cases are shown. The 
curve 2ql, reported for comparison, illustrates how the 
average current changes with the ballistic parameter I. 
In the diffusive regime (£ — > 0) the current decreases 
linearly with I (or, which is the same, with r) since scat- 
tering processes prevent carriers from gaining velocity in 
the field direction The noise in this regime corre- 

sponds to thermal (noncquilibrium) noise, and thus Sj is 
also proportional to £. Both PS schemes provide identi- 
cal results, which means that the action of the scattering 
mechanisms prevails over that of the self-consistent field 
fluctuations. When approaching the ballistic regime, the 
current increases sublinearly with I until it saturates at 
the value 0.65/ s in the ballistic limit it — > oo). Start- 
ing from I — 0.01 the noise departs from the thermal 
behavior. In the static case it obeys the classical for- 
mula 2ql typical of full shot noise (Poissonian statistics, 
no correlation among carriers is taken into account). On 
the contrary, in the dynamic case Si is strongly reduced 
with respect to 2ql as £ increases due to the correlations 
associated with the self-consistent potential. A suppres- 
sion factor of more than one order of magnitude (0.045) is 
obtained in the ballistic limit. In the intermediate quasi- 
ballistic region the dynamic noise is observed to exhibit 
a maximum value at t k 0.1. Then, the noise spectral 
density decreases with increasing £, which is somewhat 
surprising in view of a corresponding increase of the cur- 
rent. Notice that the difference between the two schemes 
starts to appear already at £ ~ 0.01, which means that 
the long-range Coulomb interactions influence the noise 
even for almost "diffusive" regimes when an electron un- 
dergoes ~ 10 2 scattering events while crossing the active 
region of the device. Another remarkable fact is that Si 
in the static case already follows the 2ql law for I = 0.03, 
although under such conditions a carrier undergoes about 
30 scattering processes in its transfer between contacts. 

To clarify the role played by different fluctuating mech- 
anisms in giving the total noise we decompose the cur- 
rent autocorrelation function into three contributions 
Ci{t) = Cv{t) + Cjv(t) + CvN(t), respectively given by 

C v (t) = ^{N) 2 (6v(t')6v(t' + t)) (la) 

C N (t) = ^{v) 2 {5N{t')5N{t' + t)) (lb) 

Cv N {t) = j^(v}{N}{5v(t')5N(t' + t) 

+ SN(t')6v(t' + t)) (lc) 

in the above equations Cy is associated with fluctua- 
tions in the mean carrier velocity Sv(t) = v(t) — (v), Cn 
with fluctuations in the carrier number and Cvn with 
their cross-correlation jyj . The spectral densities of the 
current fluctuations corresponding to these contributions 
are illustrated in Fig. 0. For £ — ► 0, i.e. in the thermal 



noise limit, Sn — > 0, Svn — * and Si — > SV, which 
means that velocity fluctuations are the main source of 
noise. Here Sy is proportional to t, thus correspond- 
ing to a diffusive behavior. It should be noted that the 
difference in the velocity-fluctuation terms Sy for the 
two schemes is practically negligible in the whole range 
of values taken by the ballistic parameter £ [Fig. 0(a)]. 
On the contrary, for the other two contributions the dif- 
ference becomes dramatic starting from I « 0.03. The 
number fluctuations increase with the ballistic character 
of transport. However such an increase is much more 
pronounced in the static case than in the dynamic case 
where, under the action of the self-consistent potential 
fluctuations, number fluctuations are significantly sup- 
pressed. The velocity-number correlations, represented 
by Svn > show an opposite behavior. Their contribution 
is positive in the static case, while it is negative in the 
dynamic case. Furthermore, for the current spectral den- 
sities calculated in the dynamic scheme, Sn and Svn 
are of opposite sign thus practically compensating each 
other. As a consequence, the current noise is considerably 
suppressed in the dynamic case with respect to the static 
case. This result reflects the fact that as carriers move 
ballistically (or quasi-ballistically) through the active re- 
gion, the dynamic fluctuations of the electric field modu- 
late the transmission through the potential minimum and 
smooth out the current fluctuations (of shot-noise type) 
imposed by the random injection at the contacts. The 
presence of scattering mechanisms, by randomizing the 
carrier velocity, reduces the fluctuations associated with 
carrier number before washing them out completely when 
the carrier motion is fully diffusive, leading the current 
noise to become thermal. 

In conclusion, within a Monte Carlo scheme we have 
investigated the influence of long-range Coulomb inter- 
action on shot-noise suppression in the transition region 
from ballistic to diffusive transport regimes. We have 
found that the strong noise suppression observed un- 
der ballistic conditions remains important also in a wide 
quasi-ballistic region, before being completely washed 
out in the diffusive regime. The role played by num- 
ber and velocity-number fluctuations in the total noise 
has been quantitatively estimated and found to explain 
the microscopic origin of the suppression mechanism. 
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FIG. 1. Current-noise spectral density Si vs the ballistic 
parameter i = Vx^rjL calculated by using static (squares, 
dashed line) and dynamic (circles, solid line) potentials. The 
dotted line represents 2qL. 



FIG. 2. Decomposition of the Si of Fig. 1 into: (a) velocity, 
(b) number, and (c) velocity-number contributions vs ballistic 
parameter I — Vthj/L for the static (squares, dashed line) and 
dynamic (circles, solid line) potentials. 
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